Shape Transitions in Hyperbolic Non-Euclidean Plates 
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We present and summarize the results of recent studies on non-Euclidean plates with imposed 
constant negative Gaussian curvature in both the Foppl - von Karman and Kirchhoff approximations. 
Motivated by experimental results we focus on annuli with a periodic profile. We show that in the 
Foppl - von Karman approximation there are only two types of global minimizers - flat and saddle 
shaped deformations with localized regions of stretching near the boundary of the annulus. We also 
show that there exists local minimizers with n- waves that have regions of stretching near their lines 
of inflection. In the Kirchhoff approximation we show that there exist exact isometric immersions 
with periodic profiles. The number of waves in these configurations is set by the condition that the 
bending energy remains finite and grows approximately exponentially with the radius of the annulus. 
For large radii, these shape are energetically favorable over saddle shapes and could explain why 
wavy shapes are selected by crocheted models of the hyperbolic plane. 
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When subjected to localized lateral swelling or growth, 
thin elastic sheets can form periodic rippling patterns 
as observed in torn plastic sheets [1], along the edges of 
leaves [2] and in swelling hydrogels Q . Since the bending 
energy of thin sheets is much weaker than stretching, the 
morphology of these structures is a result of the sheet 
bending to relieve in-plane strain. Recently, there has 
been interest in using swelling as a mechanism for shap- 
ing flat sheets into desired three-dimensional structures 

Q • This type of specified morphogenesis has been ex- 
perimentally realized by varying the radial shrinkage fac- 
tor of thermally responsive gel disks [3| and by directly 
controlling the local swelling at specific points through 
halftone gel lithography [4[. 

The equilibrium configuration of such sheets can be 
modelled as the minimum of a non-Euclidean elastic en- 
ergy that measures strains from a fixed two-dimensional 
Riemannian metric g defined on the midsurface of the 
sheet V C M 2 [5]. The intuition behind this model is 
that g encodes how the local swelling changes the intrin- 
sic distance between material coordinates and the body 
deforms to be as "close as possible" to an isometric im- 
mersion of g. In this model, deformations with vanishing 
stretching energy correspond to isometric immersions of 
g into R 3 and in the vanishing thickness limit a particu- 
lar isometric immersion, provided one exists, is selected 
by the bending energy [6|. 

This letter summarizes two studies of swelling thin 
elastic sheets that we presented in [7] and [8|. In these 
works we considered the Foppl - von Karman (FvK) and 
Kirchhoff approximations when V is an annulus with in- 
ner radius ro and outer radius R and g has correspond- 
ing Gaussian curvature Ko that is constant and negative. 
In geodesic polar coordinates (p, 6) the metric takes the 
form g = dr 2 -\-\Ko\~ 1 sinh 2 (y/|i\~o|)<i# 2 and isometric im- 
mersions of g correspond to (local) isometric immersions 
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FIG. 1: A: The evolution of the crocheted hyperbolic plane. 
As the size of the disk grows, increasing radius clockwise 
from upper left, the profile obtains a wavy pattern that be- 
comes increasingly complex. Images courtesy of Gabriele 
Meyer. B: The equilibrium shapes of hydrogel disks of con- 
stant radius R — 14mm, constant Gaussian curvature Ko = 
— .0011mm -2 , and varying thickness t — .75mm, t — .6mm, 
t — .25mm and t — .19mm clockwise from upper left [3|. 



of the hyperbolic plane H 2 [§]. 

Physical models of such sheets of have been con- 
structed by crocheting disks such that the circumference 
of the disk grows exponentially with the radius [To[- As 
the radius increases, the profile of these disks transition 
from saddle shaped to wavy with a refinement of the 
number of waves. Hydrogel disks with both radius and 
Ko fixed have also been constructed [3j]. In these gels it 
was observed that the equilibrium configurations obtain 
a periodic profile of one wavelength with a refinement of 
the number of waves n with decreasing thickness that is 
well fit by the power law n ~ t~* [3]. Motivated by these 
results we focused on deformations with a periodic profile 
with n-fold odd rotational symmetry. 

Foppl - von Karman approximation - In the FvK ap- 
proximation the deformation of the sheet is decomposed 
into in-plane e 2 x • D •— >• ^ 2 and out-of-plane displace- 
ments en : V \-> R, where e = ^/\Kq\R is a dimensionless 
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curvature. Expanding in e = ^/\Ko\R it follows that 
g ~ go + e 2 gi, where go is the Euclidean metric and 
gi = p 4 /(3R 2 ) dO 2 . Thus to lowest order the in-plane 
strain tensor is 7 = ( Vx) T + Vx + ( Vt7) t ' ~ gi > which 
retains the nonlinear contribution to the strain from the 
out-of-plane displacement. 

The energy normalized by t can then be approximated 

by 



E f [x,v}=S[ 1 ]+t 2 B[D 2 V ], 



(1) 



where r = t/(y/3Re) is a dimensionless thickness, D 2 r] 
denotes the Hessian matrix of second partial derivatives 
of 77 and the energies *S and B are the lowest order ap- 
proximate to the stretching and bending energies in [5j. 
Moreover, D 2r q is the lowest order approximation to the 
second fundamental form of the deformation and its in- 
variants tr(D 2 ri) = A 77 and det(D 2 r]) are approximations 
to the mean and Gaussian curvatures respectively. 

There are two deformations, flat and stretching free, 
that capture the scaling with r of minimizers of Ep- Flat 
deformations correspond to deformations with 77 = 0, sat- 
isfy 23 [77] = 0, and the global minimum T over such defor- 
mations can be found by solving a linear boundary value 
problem Q. Stretching free deformations correspond to 
deformations with 7 = 0, satisfy S [7] = 0, and can be 
found by solving the Monge- Ampere equation: 



det(£> 2 77) = K 



(2) 



Since Kq is constant, equation ((2]) can be solved by as- 
suming 77 is a quadratic function. This reduces the dif- 
ferential equation to an an algebraic equation. 

For stretching free deformations, Ep scales like r 2 and 
intuitively it is expected that with decreasing thickness 
minimizers should transition from flat to deformations 
that are close to being stretching free. In Q and Q we 
made this intuitive statement precise. We showed that 
the global minimum of Ep over stretching free deforma- 
tions is obtained by a saddle shaped deformation with 
77 = xy and in the vanishing thickness limit minimizers 
converge, up to rotations and translations in space, to a 
deformation with 77 = xy. 

The one parameter family of quadratic functions rj n = 
y(x — cot(7r/n)y) satisfy ((2]) as well. These functions van- 
ish along the lines = 0,7r/n and thus a deformation 
with odd n-fold rotational symmetry can be constructed 
by taking odd periodic extensions of r] n about the line 
= ir/n to form a surface 77 n . These surfaces 77 n are not 
smooth - they have jump discontinuities in their second 
derivative across their lines of inflection - but have finite 
bending energy. This is in contrast to a ridge singular- 
ity in which the bending energy diverges across the ridge 
11] . Note, that for n = 2 there is no need to make such 
extensions. 

The elastic energy of deformations with out-of-plane 
displacement 77 n scale like n 2 r 2 . Furthermore, in figure [2] 



it is clear that there is no regime in which higher period 
profiles are energetically selected over saddle shapes. In- 
deed, this statement was proved in [8] by showing that 
minimizers with n waves satisfy the following scaling for 
sufficiently small r: 
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cnr < inf £^[7, 77] < Cn r 



(3) 



where c and C are constants. This inequality captures 
the penalization paid energetically by adding more waves 
and proves that in the FvK model there can be no refine- 
ment of the number of waves with decreasing thickness. 
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FIG. 2: The normalized energy of numerical minimizers of 
Ef over deformations with n-fold odd rotational symmetry 
for r = 0.1, R = 1, K = -1 and v = 1/2. The dotted lines 
correspond to the energy of stretching free deformations with 
out-of-plane displacement 77™ while the solid line is the value 
T . The configurations plotted are colored by the approximate 
Gaussian curvature det(L> 2 77). 

Now, the exact minimizers of Ep with n waves are not 
stretching free deformations but have localized regions of 
stretching near their inner and outer radius and along the 
lines of inflection. These boundary layers are depicted in 
figure 03 Near the inner and outer radius of the disk these 
layers are the result of local stress and torque balance at 
the boundary. At the lines of inflection, the boundary 
layers are the result of corrections to the jump discon- 
tinuity in the second derivative of 77. Energetically, the 
total energy of stretching free deformations is lowered in 
these boundary layers by locally balancing the bending 
and stretching energies to 0{r 2 ). 

We showed in [8] that near the inner and outer radius 
the contributions to the bending energy arising from the 
mean and Gaussian curvatures are both relevant. Conse- 
quently, there are two overlapping boundary layers near 
the edge of the annulus in which the contributions to the 
bending energy resulting from the mean and Gaussian 
curvatures are balanced with the stretching energy. The 
intersection of these two layers form the entire boundary 
layer near the edges of the disk. The width Wk of the 
region in which the Gaussian curvature alone is reduced 
satisfies the following scaling: 

W K (0) ~ti\K \-i csc (-) cos(<9) cos (- - o) , (4) 

\ 77 / v 77 / 
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FIG. 3: The annular domain colored by the Gaussian curva- 
ture of numerical minimizers over A n for r = 0.01, ro = 0.1, 
R = 1, and Ko — —1 with n — 2,3,4,5 clockwise from the 
upper left. 



while for n > 3 the width Wh(0) of the region in which 
the mean curvature alone is reduced scales like: 



W H (0) ~ t* \K \-i^sm(^j sec(<9) sec - 6>) 



(5) 



The reason for the lack of a mean curvature boundary 
layer for the case n = 2 is a result of the surface n = xy 
being a minimal surface in the FvK approximation. This 
transition from one boundary layer in the case n = 2 
to two overlapping boundary layers for n > 3 can be 
detected in figure [3] by noticing that within this region 
K < -1 for n = 2 while K > -1 for n > 3. 

For n > 3, the width Wj of the boundary layers near 
the lines of inflection scales like 



Wl(p) ~ £3/>3 |j^ Q |- 



(6) 



In this region the bending energy is reduced by lowering 
the mean curvature by regularizing the jump disconti- 
nuity in the second derivative. This boundary layer has 
the same scaling with thickness, but not p, for minimal 
ridges formed by crumpling |ll|. Moreover, the reduc- 
tion in energy near this type of singularity is different 
from the regularization near a ridge singularity in which 
the bending energy diverges with decreasing thickness. 
Again for n = 2 there is no boundary layer since the two 
wave shape is smooth. 

Kirchhoff approximation - In the Kirchhoff model the 
deformation x : V \-> R 3 is assumed to be an exact iso- 
metric immersion of g, that is a solution of the system of 
differential equations (Vx) T • Vx = g. With this ansatz 
the energy normalized by t 3 is 



v 



AH 2 



2K 



dA 



(7) 



where H and K are the mean and Gaussian curvatures 
respectively, v the Poisson ratio of the material, and dA s 
is the area form on the surface formed by x [5, 6]. By 
Gauss's Theorema Egregium isometric immersions of g 
satisfy K = Kq and thus for isometric immersions of g it 
follows that Ek is equivalent to the energy: 

W[x] = / (k\ + k 2 2 ) \Ko\~i sinh(v^p) dpdO, (8) 
Jv 

where k\ and k<i are the principal curvatures of x. 

In this model the equilibrium configuration is selected 
by minimizing the full geometric bending energy in Q 
over exact isometric immersions. Moreover, if there 
exist isometric immersions with finite bending energy 
then minimizers of the full three dimensional energy are 
well approximated by minimizers of Ek in the vanishing 
thickness limit [6[. 

Now, isometric immersions of g correspond to solutions 
of the sine-Gordon equation: 



(j)(u,v) uv = \Ko\sm(</)(u,v)), 



(9) 



where u, v are the coordinates for a (local) parametriza- 
tion by asymptotic curves and <p is the angle between 
the asymptotic curves Furthermore, the princi- 

pal curvatures of such a parametrization satisfy k\ = 
\K \ tan 2 0/2) and k\ = \K \ cot 2 (>/2). Therefore, all of 
the geometric information is encoded in <\> and 



W[x] 
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smh(^J\Ko\p)dpde. 

(10) 

Geometrically, the coordinate curves become degener- 
ate when (j) is a multiple of tt and the principal curvatures 
diverge. At these points the surface is not different iable 
and forms sharp cusps. Therefore, isometric immersions 
of g with finite bending energy must satisfy equation (|9j) 
with the the constraint < (j) < tt. 

In we used a direct numerical scheme for minimiz- 
ing W over smooth solutions to equation (j9j). Using this 
scheme we provided numerical evidence that the princi- 
pal curvatures of smooth isometric immersions satisfy: 

m^{\h\,\k 2 \} > ^exp(2vWi?), (11) 

with the surfaces realizing this bound being geodesic 
disks lying on hyperboloids of revolution of constant 
Gaussian curvature, i.e. a saddle shape. This result 
indicates that the bending energy of smooth isometric 
immersions grows exponentially with the size of the do- 
main, a phenomenon that is not captured by the FvK 
model. Moreover, this result indicates that with increas- 
ing radius the principal curvatures of isometric immer- 
sions grows throughout the bulk of the domain. 

Furthermore, in [7[ we showed by direct construction 
that there exists exact isometric immersions A n with 
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n-wave profiles which are well approximated by rj n for 
small radii. The key to constructing these surfaces is 
the existence of a one parameter family of exact isomet- 
ric immersions Aq, called Amsler surfaces, that contain 
two asymptotic lines that intersect at an angle [l2| . 
As in the FvK approximation, if = ir/n we can take 
the odd periodic extension of the piece of the Amsler 
surface bounded between the asymptotic lines to form 
the periodic Amsler surface A n . Again, these shapes are 
not smooth but have jump discontinuities in their second 
derivatives across their lines of inflection. 



FIG. 4: The energy W of discs cut from periodic Amsler 
surfaces A n with n = 2,3,4,5 waves. The vertical dashed 
lines correspond to the radius where the principal curvatures 
of A n diverge. The dashed curve corresponds to the energy 
W of disks cut from hyperboloids of revolution. 



contradicts the known existence of smooth isometric im- 
mersions of g for finite domains [l3| . 

But, the n-periodic shapes we constructed in the FvK 
and Kirchhoff approximations do agree qualitatively with 
the experimental observations in [3]. Moreover, the 
exact isometric immersions with a periodic profile re- 
fine with increasing radius matching what is observed 
in [lOj. These observations highlight the fact that the 
system has two relevant dimensionless numbers t/R and 
e = y/\K0\R, and the models/analysis have to account 
for both of these scales, instead of just the thin limit 
t/R 0. 

Additionally, the periodic shapes observed experimen- 
tally are the result of dynamical processes. Perhaps it 
is more appropriate to model this type of swelling as a 
gradient flow. The pattern could then be selected for dy- 
namical reasons which could explain why local but not 
global minimizers appear to explain the observed shapes. 
The width of the boundary layers given by (j4][6]) are quan- 
tities that can be measured experimentally to determine 
if the local minimizers model the observed shapes. 

This work was supported by the U.S. - Israel Binational 
Science Foundation (grant No. 2004037) and the NSF 
through grant No. DMS-0807501. We would also like to 
thank Gabrielle Meyer for providing photographs of the 
crocheted model of the hyperbolic plane. 



Isometric immersions of geodesic balls can be "cut" 
out of A n . In figure |4] we plot the energy W of such 
immersions and compare their energy to those cut from 
hyperboloids of revolution. To be precise, there is a one 
parameter family of hyperboloids of revolution and for 
each radius we selected the hyperboloid with the lowest 
value of W. This figure illustrates that for large values 
of e = ^/\Ko\R the n-wave shapes are energetically pre- 
ferred over their smooth counterparts. Additionally, in 
[7[ we showed that there is a critical radius R n ~ ln(n) 
where the principal curvatures of A n diverge. Beyond 
this radius A n no longer has finite bending energy and a 
profile with a higher number of waves is selected. This 
gives a natural mechanism for the refinement of the num- 
ber of waves with increasing radius as observed in [lOj . 

Discussion of results - From the studies in [7[ and [8[ it 
is clear that the refinement of the number of waves with 
decreasing thickness observed for the gels in [3| cannot 
be directly explained by global minimizers in the FvK 
and Kirchhoff approximations. Specifically, we proved in 
and @ that for these gels there are only two types of 
global minimizers - saddles and flat sheets - in the limits 
of thin and thick sheets, and no intermediate asymptotic 
regime in which shapes with more waves have lower en- 
ergy. Moreover, the bending content of the configura- 
tions in [3| diverges according to the scaling t~ x which 
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